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We demonstrate that the quantized U(l) lattice gauge theory in four Euclidean dimensions can 
be obtained as the long time average of the corresponding classical U(l) gauge theory in 4+1 
dimensions. The Planck constant Ti is related to the excitation energy and the lattice constant of 
this classical template. 
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I. INTRODUCTION 

The nonlinear classical dynamics of gauge fields is 
known to be strongly chaotic [HQ]- The classical SU(2) 
gauge theory defined on a three-dimensional lattice has 
been studied numerically in considerable detail and was 
shown to be a globally hyperbolic (Anosov) system 
For such systems, any generic initial gauge field config- 
uration wanders ergodically in the phase space of field 
configurations. This motion leads in the infrared limit 
to a stationary distribution of lower dimensional config- 
urations [4(. These, as a result of the higher dimensional 
classical dynamics, are distributed exactly as required by 
the vacuum state of the Euclidean quantum theory on the 
lower dimensional space. 

This phenomenon, which has been called chaotic quan- 
tization can be simply considered as a physical 
realization of the method of stochastic quantization [8| . 
It is a consequence of two deep relationships: (1) that 
the long time properties of certain deterministic classical 
systems can be described by the methods of statistical 
physics 0, and (2) the correspondence between classical 
statistics and Euclidean quantum mechanics. In the ap- 
pendix we demonstrate this correspondence for the quan- 
tum harmonic oscillator in detail. 

According to the chaotic quantization scenario it is 
necessary for the higher dimensional (second time) dy- 
namics to be chaotic in order to appear as a quantum 
theory in the lower (one-time) dimensions. From the in- 
vestigation of the U(l) gauge theory presented in this 
article we shall see that the existence of the lattice is 
essential for this scenario to work in those cases when 
the dynamics would not be chaotic in the weak coupling 
limit. The chaotic quantization may fail for systems not 
being chaotic on finite lattices as well as in case of ob- 
serving too short second-time of the dynamical evolution. 
This opens up towards the possibility that not all inter- 
actions, in particular gravity, have to be quantized in 
nature. 

This way the chaotic quantization of gauge theories 
also lends support to 't Hooft's proposal 7] that quan- 
tum mechanics could arise from an underlying dissipa- 
tive, yet microscopically deterministic classical dynam- 
ics. We will not explore this speculative avenue here and 
confine ourselves to reporting our numerical evidence for 



a correspondence between the classical U(l) gauge theory 
in 4+1 dimensions and the quantum U(l) gauge theory 
in 4 dimensions. 

The classical compact U(l) gauge theory without mat- 
ter fields is the simplest gauge theory still exhibiting 
chaotic dynamics in calculations on three-dimensional 
lattices ||- Although its continuum limit is not chaotic, 
in contrast to non-Abelian gauge theories where chaos 
survives even in this limit, the classical dynamics of the 
compact U(l) lattice gauge field is chaotic at finite lattice 
spacing a and sufficiently strong coupling g (ga 2 finite) . 
As in the case of non-abelian gauge theories few mode 
models, like the xy-model ^2| can be studied, a 3-link 
section of U(l) lattice gauge theory reveals a system of 
harmonically coupled pendulums, easy to see numerically 
drifting in chaotic behavior at high energy. The dimen- 
sionless parameter controlling the strength of chaoticity 
of the lattice gauge theory is g 2 aE p , where E p is the av- 
erage total energy per elementary lattice plaquette. 

Our conjecture of the correspondence between the 
(4+l)-dimensional classical gauge theory with chaotic 
dynamics and the 4-dimensional Euclidean quantum field 
theory can be expressed in the following relation between 
the 4-dimensional Planck constant h and two physical pa- 
rameters of the higher dimensional classical theory: the 
temperature T and the lattice spacing £|: 



h = aT. 



(1) 



Our philosophy here views h as a constant of nature fac- 
torized into two underlying properties of the true (4+1)- 
dimensional world. An analogous case is the relation of 
classical electrodynamics factorizing the speed of light c 
into electric and magnetic properties of the vacuum: 



1/c 2 = eoMo- 



(2) 



Taking the invariance of c as an axiom in the theory of 
special relativity, one can derive the Lorentz transforma- 
tion law without the need for any reference to electric or 
magnetic fields. Maxwell theory on the other hand, as 
a classical field theory, regards the dielectric constant eo 
and the magnetic permeability /j,q as independent prop- 
erties of the physical vacuum. Light waves are solutions 
of the Maxwell theory, and the speed of light is calcu- 
lable. The proposed relationship between quantum field 
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theory and an underlying classical field theory in a higher 
dimension is analogous. 

In this letter we report the results of numerical sim- 
ulations of a U(l) lattice gauge theory both in 4 and 
4+1 dimensions. The 4-dimcnsional Euclidean quantum 
gauge theory was simulated by standard quantum Monte- 
Carlo techniques. The 4+1-dimensional classical theory 
was simulated by numerical solution of the differential 
equations describing its Hamiltonian dynamics. In both 
cases, 4 4 hypercubic lattices were employed. Lattice links 
beginning at point x and pointing in the /x direction arc 
associated with phases, A^x), of unimodular complex 
numbers U x ^ — ex^>[igaA^{x)] € U(l) in these theories. 



II. CLASSICAL AND QUANTUM LATTICE 
GAUGE MODELS 

In order to discuss the lattice regularization of classical 
and quantum field theory in parallel, we use a common 
notation as far as possible. Our starting point is the def- 
inition of the link variables U , which are always elements 
of the local U(l) group, but the interpretation, and hence 
the decomposition, of the phase can be different. In or- 
der to relate a physical interpretation to the lattice model 
phase, gaA, we consider the product of charge and vec- 
tor potential, qA (in the CGS system with c = 1), which 
is the interaction energy associated with the field living 
on a link of length a with a physical test charge q. This 
expression must be divided by another quantity with the 
dimension of an energy to generate a dimensionless phase 
variable. This quantity, the unit link energy e, is used as 
an energy "standard" , but it does not appear in the equa- 
tions of motion of the theory. Also, we are free to choose 
e differently in the quantum and in the classical theory. 
The general link variable is now given by 

U L =exp(iqA/e). (3) 

In the standard quantum lattice gauge field theory the 
energy scale associated with a link of length a is 

e Q = h/a. (4) 

In a classical theory, however, no reference to the Planck 
constant h is allowed. In this case we will use instead the 
Coulomb energy associated with two charges <zk at the 
endpoints of the link of length a, which in three spatial 
dimensions is given by 

eK = <?K/a- (5) 

By comparing the classical (K) and quantum (Q) sys- 
tems of U(l) lattice gauge fields, we shall consider con- 
figurations characterized by the same link phase gaA. 
Using a common lattice spacing a and a common cou- 
pling strength g, this convention naturally leads to the 
consideration of the same vector potential field, A, in the 



continuum limit. We can express this correspondence by 
the relation 

<?Q/ e Q = W £ k = ga. (6) 

As a consequence the classical charges as reference 
sources for the unit link energy are interpreted differ- 
ently in the two cases: qQ = gh in the quantum the- 
ory and <?k = 1/5 in the classical lattice theory. In a 
sense, the two approaches are dual to each other, with 
QQQk = Ti, i.e. a high value of the charge corresponds to 
strong coupling in the quantum, but to weak coupling in 
the classical theory. 

We note that the conventions chosen by us here are 
not unique. Alternative conventions, e.g. relating ck to a 
classical charge in four spatial dimensions, as 6k — q 2 /a 2 
are also possible, but they would not lead to different 
conclusions. In the end, the classical theory as well 
as the quantum theory are defined through the interac- 
tions of the dimensionless link variables U. The classical 
Hamiltonian Hk[U], being scale invariant, is a nonlinear 
function of the link variables, multiplied by a constant 
parameter carrying the dimension of energy. Likewise, 
the action of the quantum theory Sq[U], which depends 
on the scale only through the gauge coupling l/g 2 (a) as 
an overall factor, has the dimension of an action. If, as 
we will show below, the two theories are related to each 
other, this implies a relationship between the dimension- 
ful factors of Hk and Sq, which does not depend on the 
conventions used in the definition of the phase of U. 

Physical quantities are related to the oriented prod- 
uct of link variables around an elementary square, a pla- 
quette. Using the notion of lattice forward derivatives, 
ad^f = f(x + flEp) — f(x), (where e M is a unit vector 
in the corresponding direction) the plaquette phase sums 
are associated with the local rotations of the vector po- 
tential, with the field strength tensor 

F„ v {x) = d^A v {x) - dvA^x). (7) 

An elementary plaquette variable is related to a compo- 
nent of the field strength tensor: 

U x ^ v = cxp [iga 2 F^x)) , (8) 

and the lattice sum over all plaqucttcs 

E P = ££(l-Retr(t/ x ,^)) (9) 

is related to the physical energy or action, depending on 
the dimensionality of the considered space. For the U(l) 
group, 

1 - Re tr(U x ^) = 1 - cos (ga 2 F^{x)) , (10) 

which in the continuum limit a — > approaches 
(l/2)g 2 a 4 (F^) 2 . The plaquette sum approximates 

Sp^V^^x). (11) 
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Over a 4-dimensional lattice it is proportional to the ac- 
tion, Sp = g 2 Sii and can be used to express the loga- 
rithmic weight of a configuration: 



Si/h = £ 
9 ri 



(12) 



where (3 is a dimensionless constant playing the role of 
a (fictitious) temperature. Note that the fine structure 
constant can be expressed as a = q 2 /4irh, yielding the 
familiar relation 



g 2 h 2 
Ann 



Air 



1 

4^/3 



TZZ (13) 



for quantum lattice gauge theory, and 



1 



or 



Airhg 2 



4tt' 



(14) 



for the classical theory. The duality of these interpreta- 
tions is manifest here, as well. 

In the (4+l)-dimcnsional classical lattice gauge the- 
ory the summation over all plaquettes of a 4-dimensional 
spatial lattice is proportional to the (four-dimensional) 
magnetic energy. Since the classical interpretation can- 
not refer to the Planck constant, we arrive at the follow- 
ing expression for the energy by using ek as the energy 
unit: 



1 

9 a 



(15) 



Here we used (JBJ to eliminate ok from the definition (JSJ) of 
ex- In the continuum limit, the magnetic energy becomes 



,4 V" 



1 



^ 4o 

X,fXV 



(16) 



confirming that the physical gauge fields in 4+1 dimen- 
sions, 

,4(4+1) 

arc related to those in 4 Euclidean dimen- 
sions, by a factor a -1 / 2 . 

We now come to our main point. The conjecture 
h = aT is supported if the quantum and classical lat- 
tice simulations can be brought into a correspondence 
relating configuration weights, and thereby all physical 
expectation values, by 

exp{-S 4 /K) oc exp(-^ nagn /T) oc exp(-/?£ P ). (17) 

We have used proportionality signs in this relation, be- 
cause the overall normalization of the weights is not rele- 
vant. Also, the relation l|17|l defines the temperature T as 
the equipartitioning energy for the magnetic energy of the 
classical lattice field. This designation obviously requires 
that the classical theory is ergodic, so that a single field 
trajectory generates configurations with a microcanoni- 
cal equilibrium distribution. 

Quantum Monte Carlo algorithms generate configu- 
rations of lattice link elements according to the Gibbs 



weight w[U] oc e~ /3Sp = e~ s \- u ^ h . The classical Hamilto- 
nian approach, on the other hand, divides the Hamilto- 
nian into its electric and magnetic parts: 
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^ 2a 

X . fJi 



H = a' l y ' —A z + E masn [U}. 



(18) 



This formula leads to the following dimensionless Hamil- 
tonian lattice model, where the time is measured in lat- 
tice spacing units a and time derivatives take the dimen- 
sionless form a(d/dt): 



g 2 aH = i tr (aUallA + S P 



(19) 



Now the link variables U x „ are still defined on a four- 
dimensional lattice, but they are functions of an addi- 
tional, continuous and scaled time variable t/a. The evo- 
lution of the lattice configuration occurs in this fifth di- 
mension, usually called "fictitious" time in the context of 
stochastic quantization. In the present context, however, 
the time-like dimension is considered as physical, albeit 
unobservable at long distances and low frequencies. One 
reflection of this difference is that we do not add an exter- 
nal heat bath or white noise to the classical dynamical 
equations; we here consider pure classical Hamiltonian 
dynamics. 



III. TECHNICAL ASPECTS OF THE 
SIMULATIONS 

We now explain some technical aspects of our simula- 
tions, before we report the results. For the purpose of 
generating generic initial field configurations, and as a 
matter of convenience, we prepared the initial configura- 
tion for the classical Hamiltonian simulation by Monte- 
Carlo "heating" on a five-dimensional lattice. The /?5 
parameter chosen for this simulation determines the av- 
erage value of the action S$ and of the five-dimensional 
electric and magnetic energies. The selected configura- 
tion was then converted into the initial data for a four- 
dimensional, real-time lattice calculation as follows: The 
space-space-like plaquettes and links were taken with- 
out change from the four-dimensional hypercube located 
at X5 = 0. The periodic boundary condition in the 
xs-direction allowed us to identify the link phase for 
x^/as = — 1 with those at 3:5/05 = N$ — 1 for the con- 
struction of the initial values for the time derivatives U 
on the hypercube in terms of the link variables attached 
to the plane x§ — 0. This construction also used the 
link phases in the x§ direction to make the time deriva- 
tives gauge covariant. By construction, the U variables 
then are orthogonal to the corresponding U variables, as 
required by the Hamiltonian time evolution. We now 
present the details of our algorithm for the construction 
of the initial gauge field configurations. 

We denote by U+ the triple product of link variables in- 
cluding the value at 2:5/05 = 1 (the first argument refers 
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to the site coordinates x = {x\, ■ ■ ■ , £4}; the second ar- 
gument indicates the value of 2:5/05): 

U+ = U 5 (x, 0)U x (x, l)Ul(x + a s , 0) (20) 

and similarly by U- that product including the previous- 
time value 

U- = U$(x, -l)U x (x, -l)U 5 (x + a s , -1) (21) 

Both triples close a plaquette with the U = U x (x, 0) link. 
These plaquettes encode the electric field values at x$ ± 
as/2. Using the scalar product notation for U(l) group 
elements: 

(A,B) = cosO/u - (22) 

we can determine the electric field energy densities 
s± = ^E x (x, ±^f) as follows: 

(U+,U) = l-(ar a s ) 2 e+, 

(U-,U) = l-(a 5 a s ) 2 e_. (23) 

We also define an interpolating electric energy density e 
by using a double-sized plaquette: 

(U+,U-) = l-{2ar a s fe. (24) 

The momentum P stands for Udt in the classical cal- 
culation. We relate U to the link variables of the five- 
dimensional lattice by identifying dt — 05 for the initial 
state of the Hamiltonian simulation. Later on we use 
this initial value of P and U for updates in much smaller 
steps than the original lattice spacing, dt/a 5 <§; 1. (Typ- 
ical values are 0.01 and 0.001.) The initial P can be 
expressed as a linear combination of the U and U± link 
variables; to leading order it is a difference (with a pos- 
sible admixture of second derivative), 

P = \{U + -U-) + 5 -{U + + U--2U). (25) 

The parameter 8 is obtained from the orthogonality con- 
straint 

1 8 

(P, U) = -(a 5 a s f(-e + + e_) - -(a 5 a s ) 2 (e+ + e_) = 0. 

(26) 

This equation can be easily solved for 8, and P is obtained 

as 

p = e-U+-s+U- - (e_ -e+)U 
e-+e+ 

The electric energy per link is related to the dimension- 
less kinetic term 

\(P,P) - (28) 

The lattice simulation programs were coded in C + + 
constructing classes for arbitrary sized and arbitrary di- 
mensional lattices with an established site, link, plaque- 
tte and elementary cube index system. In particular in- 
dices for calculating gradient, rotation and divergence 



were taken care of. Classes for group variables were also 
constructed with group multiplication, Haar measure and 
random generation functions. In case of U(l) a single real 
phase represented the group variable. Data i/o handling 
was coded in compressed binary form with a standardized 
header of lattice size information. 

Algorithms for cold, hot and exp(— (3(1 — cos ip)) dis- 
tributed initialization, as well as Monte Carlo heating 
algorithms based on the original Metropolis rejection 
method or in another version on a local heat-bath mech- 
anism were implemented for producing lattice configu- 
rations. Functions like the complex Polyakov line aver- 
age over the lattice space volume and the action were 
included. A conversion program between data on a d 
dimensional lattice and a Hamiltonian pendant of the 
configuration on a (d — 1) dimensional lattice (plus an 
almost continuous time variable) determined U and U 
values to begin a classical solution of the equations of mo- 
tion. These were solved by a program obtaining updates 
in a stabilized half-time-step algorithm, ([/<—£/ + P/2, 
P <- P + (dt 2 )F, U <— U + P/2 again), which ensures 
that energy is conserved up to a precision of 0(dt 3 ). We 
have used typically dt/a = 0.01. 

We now give a brief summary of our tests of these pro- 
grams for performance. The group multiplication, trace, 
determinant and Haar-measure routines were tested on 
special (unity, exp(i7r)) and on random phase elements. 
The random generators, one using the simple rejection 
method with 1 as majoring function and one using a 
Gaussian to decide whether to reject or not, were tested 
for averages, variances, auto-correlation and computa- 
tional time. The transient effects by Monte Carlo heating 
were monitored by the average action, by the Polyakov 
value squared and by the magnetic monopole density. An 
optional acceleration of heating was implemented by the 
link phase mirroring technique inserted between the tra- 
ditional heat-bath steps. The usual phase transition as a 
function of the quantum Monte Carlo coupling constant 
has been performed at the expected value for 4 4 , 6 4 , 8 4 
and 12 3 x 4 lattices. In the classical equation of mo- 
tion solver the energy conservation and the equipartition 
of the kinetic and potential energy parts has been moni- 
tored. All averages presented in the following are already 
screened for transient effects, which have been removed 
by eliminating early parts of the trajectory. 



IV. NUMERICAL RESULTS 

Microcanonical equilibration of classical Hamiltonian 
lattice fields has been studied earlier for the non-Abelian 
groups SU(2) and SU(3) in three spatial dimensions y]. 
Here we repeat this study for the U(l) system in Ad in or- 
der to make sure that this system also approaches a ther- 
mal distribution of plaquette energies during its chaotic 
evolution. Fig^ plots the frequency of occurrence of the 
r = 1 — cos(ag 2 F) plaquette values in the allowed inter- 
val of [0, 2] on a logarithmic scale for a corner plaquette 
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at different times (sampied at each 5-th time-step of an 
evolution of 100.000 steps of each having dt/a = 0.01). 
Due to the periodic boundary conditions there is noth- 
ing special about this position. One observes that the 
distribution of plaquette values is well approximated by 
an exponential law multiplied by the inverse square root 
function (red histogram). We note by passing that a 
similar distribution has been observed in earlier studies 
of 3-dimensional SU(3) lattice g aug e models, but with a 
prefactor of r 3 instead of r -1 / 2 flOj | . 

For the sake of comparison the same distribution is 
plotted for the 4-dimensional quantum Monte-Carlo algo- 
rithm at inverse coupling f3 — 1.3 (blue histogram). The 
correspondence is made due to similar values of the aver- 
age Polyakov line squared {{P 2 )mc ~ (-P 2 ) chaos ~ 0.44) 
at close values of average plaquette traces (g 2 Si ~ 0.225 
and g 2 aE™ a9n rj 0.23). The agreement between these 
two distributions, not only in the shape, but also in abso- 
lute number, already demonstrates that we are producing 
instances of lattice configurations with the same proba- 
bility by applying either the traditional quantum Monte 
Carlo or by the chaotic dynamics method. 




O 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

r = 1 - cos(ag 2 F) 



FIG. 1: The distribution of plaquette values is well ap- 
proximated by an inverse square root times an exponential 
decay law. This provides numerical evidence for the self- 
thermalization of the classical U(l) lattice system due to its 
chaotic dynamics (red histogram), as well as of the Monte- 
Carlo heat bath algorithm (blue histogram). 

We now further demonstrate the validity of the relation 
(|17fl by showing the equivalence of the two theories for an 
important observable, the expectation value of the modu- 
lus squared of the Polyakov line on the respective periodic 
lattices. We note that this is the standard order param- 
eter of the lattice gauge theory, which vanishes at strong 
coupling and is nonzero below a critical coupling /3 C w 1, 



indicating the transition to the Coulomb phase. Figure[3] 
shows the Polyakov line modulus squared averaged over 
the 3- volume of the 4-dimensional lattice and over 20,000 
quantum Monte-Carlo configurations (red squares) as a 
function of the average 4-dimcnsional lattice action per 
plaquette £p = g 2 Si. 

On the same plot the same order parameter is shown 
(green dots) as a function of the expectation value of the 
lattice magnetic energy Sp = g 2 aE m&g , per plaquette, of 
the classical configuration. The average here was calcu- 
lated, after an additional brief Hamiltonian equilibration 
on the classical four-dimensional lattice, as an ergodic 
average by temporal sampling of a single evolving lattice 
field configuration. We note that these two sets of points 
follow the same scaling law. 

The reason that we do not plot the modulus of the 
Polyakov line as a function of the inverse coupling /?, 
as it is usually done, is that the classical Hamiltonian 
(|T%|l does not contain such a coupling constant. Its so- 
lutions are solely characterized by the value of the total 
energy or, because of its ergodicity properties, by the av- 
erage value of the magnetic energy. The fact that our 
simulation points obtained for the quantum field the- 
ory on the four-dimensional Euclidean lattice and the 
(4+l)-dimensional classical theory coincide when we im- 
pose the relation H17|) constitutes the desired evidence. 
(One might consider relating the two results via the cou- 
pling constant (3§ used in the generation of the initial 
conditions for the classical field configuration. However, 
this would be inappropriate, because the five-dimensional 
lattice provides us solely with a convenient algorithm 
for generating randomized initial data. We could have 
started the classical calculation with an arbitrary initial 
configuration of the same total energy, without any ref- 
erence to a five-dimensional lattice, and would have ob- 
tained the same results, albeit after much longer micro- 
canonical equilibration.) 

In order to further confirm the validity of the rela- 
tion between the four-dimensional quantum and (4+1)- 
dimensional classical lattice U(l) theories it is illustrative 
to look at scatter plots of sampled values of the Polyakov 
line in the complex plane. In Figure |21 the results from 
the four-dimensional quantum Monte-Carlo simulation 
are shown as red dots, and the results from the (4+1)- 
dimensional classical Hamiltonian evolution for a single 
trajectory are represented by green dots. The correspon- 
dence in the four parts of Fig. |3| is again made by se- 
lecting pairs of simulations satisfying the relation l|17|l . 
The overlap of the distributions is excellent, both for the 
modulus and the phase of the Polyakov line. Some initial 
points in the center of the rings for supercritical couplings 
/?4 = 1 / <?4 > 1 are artefacts from an initial heating phase 
in the quantum Monte Carlo simulations. They would 
be absent, if we had started the sampling process at a 
later time. In the classical Hamiltonian evolution the er- 
godic sampling was started also only after a transitory 
period, in which the electric and magnetic field energies 
equilibrated (in 4+1 dimensions, the ratio of electric and 
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\0 5-dim eom 



CM 



E 5 / T = S 4 / h 



4-dim MC 



0.2 0.4 0.6 0.8 1 

g%E 5 ™°", g 2 S4 

FIG. 2: The order parameter, the absolute value square of the 
Polyakov line averaged over the lattice and over many config- 
urations obtained by classical Hamiltonian dynamics is plot- 
ted against the magnetic energy of the four-dimensional spa- 
tial lattice (open circles), and the same quantity obtained by 
traditional quantum Monte-Carlo simulation on a Euclidean 
four-dimensional lattice against the four-dimensional action 
(filled boxes), respectively. These results coincide, if and only 
if E 5 /T = S 4 /h. 



magnetic energy is 2:3, not 1:1 as in the usual (3+1) 
dimensional case). 



basis of non-Abelian gauge theories in 3 and 3+1 dimen- 
sions - also works for the compact Abelian lattice gauge 
theory in 4 and 4+1 dimensions. The correspondence be- 
tween the classical and quantum gauge theories is given 
by the general formula Ti = aT, which encodes physi- 
cal properties of the higher dimensional theory into the 
Planck constant. Our result supports the speculation 
that chaotic quantization may be a physical mechanism 
for the quantization of (gauge) fields. This could make 
it eventually possible to construct a classical field theory 
encompassing both, gravity and the standard model of 
particle physics. In this framework, Planck's constant h 
would become a parameter of the low-frequency, long- 
second-time limit of the fundamental classical field the- 
ory. 



Finally we would like to address the question whether 
factorizing the Planck constant would be tantamount to 
the construction of a hidden variable theory. We be- 
lieve that this is not the case, since none of the estab- 
lished rules of quantum physics are violated. The four- 
dimensional quantum field theory is given by a Euclidean 
functional integral, which is exactly the path integral that 
defines the vacuum sector of the compact U(l) gauge the- 
ory in (3+1) dimensions. The higher dimensional clas- 
sical dynamics "acts" as a quantum field theory in four 
Euclidean dimensions. 



V. CONCLUSION 

In conclusion, we have demonstrated that the mecha- 
nism of chaotic quantization - conjectured earlier on the 
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APPENDIX A: CHAOTIC QUANTIZATION OF 
HARMONIC OSCILLATOR 

In order to help to gain a deeper insight into the chaotic 
quantization mechanism we present here a limiting case, 
the harmonic oscillator. Of course, pure lattice gauge 
theories cannot be used as models of the oscillator be- 
cause they lead to massless particles in the continuum 
limit with a trivial dispersion relation u> — in the + 1 
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dimensional case of ordinary quantum mechanics. An ad- 
ditional ingredient is needed for this purpose; we choose 
a complex scalar (Higgs) field which generates mass for 
the U(l) fields leading to the desired harmonic oscillator 
description in the low dimensional case. We consider the 
following action, 



S = 



J (-jF^ + l\D^\ 2 -vm) 



d A x 



(Al) 



where the Higgs potential V prefers a non-trivial vacuum 
state |$(t, x)\ — R. The Euclidean lattice version of this 
theory is investigated on a narrow, N x 1 stripe in the 
t — x plane (here t denotes the "ordinary" Euclidean time 
direction), having a periodic boundary condition in the 
dummy x direction. The link phases starting at sites 
tj = ja and pointing in the (extremely short) x direction 
are denoted by ipj, the ones pointing in the t direction 
by pj. The complex fields at site positions are given by 
<&(tj) = Rid 1 i. The lattice action is decomposed as 



with 



SE.lat — Sg + S x + St, 



(A2) 



S x 

Si 



2 ^ 3 

j 



a 2 \R j+ xe i ^+ 1+p ^ - Rje* 9 * 



(A3) 



In the Higgs phase with a Lorentz-invariant (in the case of 
+ 1 dimensional quantum mechanics, just simply static) 
modulus Rj = R, terms containing tpj decouple from the 
rest. This constitutes the lattice action for the harmonic 
oscillator, reducing the dynamics to that of a long chain 
of minimal plaquettes in the t direction: 

SiatW\ = ^^(l-cos(( y 5j+i-(p i ))+(a-R) 2 5I( 1_cos ^)- 

9 j 3 

(A4) 

With the correspondence u> = gR and ipj = gqj yjrnja we 
get the familiar action of the harmonic oscillator in the 
small phase (either weak coupling, or small amplitude, 
or small oscillator mass) limit 



Siatk) = a 22 




(A5) 



For a — > this leads to the continuum Euclidean action 
of a single harmonic oscillator: 



*w - / 



q 2 dt. 



(A6) 



We now turn to the analysis of the Hamiltonian descrip- 
tion of this lattice chain model l|A4|l describing its classi- 
cal evolution in a second time s. We use the Hamiltonian 



— Slat[ 

a 



(A7) 



The classical equation of motion, presented in a scaled 
time variable s/a becomes 



V., 



d(s/a,y 



sm(</5,-+i -tpj)- sm{ipj - (pj-i)- (au) 2 



in ipj. 

The only parameter in this s-dynamics is the scaled fre- 
quency G = au. In a molecular field approximation each 
link phase (or plaquette phase) would behave like a pen- 
dulum with a harmonically perturbed pivot point. Such 
systems are known to be chaotic [Tlj . 

In order to understand how the classical chaotic me- 
chanics leads to a distribution in the second time s, 
which resembles that distribution applied in Euclidean 
quantum field theory, the key property is the relation 
between the equations of motion and a Fokker-Planck- 
Kolmogorov equation. This relation has also been stud- 
ied at length in chaotic dynamics and the circumstances 
when the chaotic motion can be replaced by a diffusion in 
phase space governed by a Fokker-Planck type equation 
have been investigated. Here we follow a presentation 
given in Chapter 6 of Ref . [lll| . 

First the differential equations are solved in discrete 
time steps, As, generating a mapping in the phase space: 



pi 



Fj As 



(A9) 



with the following force in the lattice chain model 



F; 



m d 



a d(p : 



Slo 



(A10) 



The mapping should be in action-angle variables; the lat- 
tice model is readily formulated in terms of the angles 
ipj. The corresponding scaled action variable we denote 
by Ij = a(ckpj/ds). The (2-time) Hamiltonian contains 
the following kinetic energy term 



Hi 



= m v ( —Y = 1 1 v j 2 

2 ^\ds ) 3 2 a2^ ! 



(All) 



allowing us to write the mapping in dimcnsionless action 
- angle variables as follows 



V j = tpj + 



-I 



As 

a 
A.s 



g q Slat- 

a oifij 



(A12) 
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An equivalent Fokker-Planck-Kolmogorov equation is 
then valid for the distribution of 7j in a long-term sam- 
pling of phase space trajectory: 

£--&**> + 55555 (Bb ' 5 (A13) 

with the following drift and diffusion coefficients 
Ai = ^((Ji -/<)), 

Bij = (A14) 

Here the averaging over different observation moments s 
and different trajectories (wherever started) can be re- 
placed by an average over the angle variables. This is 
a central statement of chaotic dynamics, giving a statis- 
tical physics perspective to deterministic dynamical sys- 
tems. This is as well the key mechanism beyond the 
chaotic quantization, arriving at a Fokker-Planck corre- 
spondence, which used to be the starting point of the 
stochastic quantization method. 

For our lattice oscillator model we get A, = and 

% = v ( (1 + HT )( ^ " \ {5 ^ + 5i+1 ^) ■ 

(A15) 

meaning that the diffusion matrix, the coefficient of the 
scaled time-step As /a in the matrix is readily recog- 
nized as the the harmonic oscillator matrix in the quan- 
tum mechanical Euclidean path integral formalism. 

At this point we still seek to determine the long second- 
time limit of the diffusion matrix B^ . This can be under- 
stood with the help of the Green-Kubo formula j(| which 
connects the microscopical dynamics with the linear re- 
sponse approximation. The diffusion-like Fokker-Planck 
equation is actually the leading order cumulant expan- 
sion to the more general case, the Fourier transform of 
the probability is taken in the infrared limit associated 
to the basic action variables: 

T k = (e lk -^) a exp (-fc 4 {A/ 4 AIj)kj) . (A16) 

Back Fourier transformation leads then to the second 
derivative according to the Ij variables; the first deriva- 
tive (drift) term vanishes for time reversal dynamics, such 
as governed by conservative Hamiltonians. The diffusion 
coefficient in scaled time s/a becomes 

Bij = («^7 A ^) = (A17) 

Assuming now that the correlation of the action vari- 
ables depends only on the absolut value of second-time 
difference, as it is the case both in stochastic and chaotic 
processes, we arrive at 



The characteristic correlation of the force (I) is an ex- 
ponential decay both in the stochastic, as well as in the 
chaotic quantization. In the former case a Langevin equa- 
tion with white noise leads to forgetting, in the latter case 
the instability of nearby trajectories expressed by posi- 
tive Lyapunov exponents. The general behavior is then 

(/<(*') j 3 -(0)> = e -^'(/ i (0)7 J -(0)). (A19) 

The second-time integral in IjAf 8(1 can now be analyti- 
cally done leaving us with 

B^ = (A20) 

For a short s-time evolution it gives back the result of 
(|A15fl . in the long term, however, leads to a constant 
times the matrix describing the oscillator's path integral 
in ordinary quantum mechanics: 

km B id = — ((1+ ^-)S ij - hSt-u + S i+1J ) \ . 
s-^oo aj \ Z Z J 

(A21) 

This proves our conjecture for the harmonic oscillator, 
without extended numerical simulations. 

Finally some remarks are in order to explain how di- 
mensionful scales and parameters enter in the interpreta- 
tion of lattice model results. Of course, both the familiar 
quantum Monte Carlo and the second-time Hamiltonian 
equation of motion (shortly EOM) method deals with 
scaled quantities on the lattice (and with a scaled time) 
not having any length (fm) or energy (MeV) dimensions. 
Besides the lattice spacing a, giving the unit of length, 
either the unit of quantum action, h a 197 MeVfm is 
given, or in the EOM method instead of ft an energy, 
which is conserved in total, is given. 

The study of the diffusion coefficient in phase space 
describes the chaotic evolution in mid terms, the elapsed 
second time As to be taken between an autocorrelation 
time s c when initial phase points are remembered by the 
evolution trajectory (there is an exponential forgetting 
of nearbyness of the initial conditions set by the leading 
Lyapunov exponent) and the longer time when equipar- 
tition completes. As in the case of simple Brownian mo- 
tion, the kinetic energy starts to grow linearly for short 
time s but levels off at an equipartition value, in our case 
we have 

<2^£ ? If) = N "2 {1 - 6 ~ 7S) ~* NT/2 > (A22) 

where N is the ordinary time extension of the path inte- 
gral lattice. (This is the correct scaling of Brownian mo- 
tion parameters with "volume".) Expressing the prod- 
uct of the equipartition temperature T and the inverse 
equipartition time 7 in a lattice scaled fashion we arrive 
at 



B^ = I ' (ii(s')ij(0))ds'. (A18) 
Jo 



(A23) 
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containing the average eigenvalue of the diffusion matrix. 
On the other hand the long time equipartition reaches 

(\s lat ) - \TN, (A24) 

while the distribution T approaches the canonical one 

T oc e~ N = e -^/ T = e -2<W(«T)> (A25) 

This canonical distribution is eventually interpreted as 
the Euclidean quantum wave function, T ~ exp(— S/K), 
(continuing back to real time it would be exp(iS/K) con- 
stituting the solution of the Schrodinger equation). This 
interpretation of the chaotic distribution allows us to 
connect the Planck constant experienced in the 1-time 
quantum mechanics with parameters of the phase space 
diffusion of the 2-time Hamiltonian: 

n=\aT=^XB). (A26) 

For the oscillator in the N — > oo limit the average eigen- 
value becomes X(B) = (2 + (aui) 2 ) leading to 

1 1 / (auj) 2 \ /A „„, 

9 2 aj V 2 / 
I 



or with a — > at finite g and = gR values, the same 
value as for free fields, 



1 f 

g 2 -fa' 



relating Planck's constant with the leading Lyapunov ex- 
ponent of the chaotic dynamics in s-time. 

For the classical simulation input is the scaled energy, 
g 2 aE 5 which is conserved, and the lattice size N. For 
the quantum MC simulation input is (3 = \/(g 2 Ti) and 
the size N . From the above equipartition by chaotic dy- 
namics we got g 2 U — 1/(07) — 9 1 'aE 5 / '(2N) . Since by 
construction the half of E 5 is S\/a at equipartition, we 
tend to a final state where (Si at ) = Nh . As we equiparti- 
tion energy in 5d, making a temperature T, with proper 
scaling we create an action in Ad with an Ti value per 
degree of freedom. In this sense the value of h is deter- 
mined by the scaled total energy per degree of freedom 
into the classical EOM simulation. Its interpretation and 
descendence are, however, quite untraditional in chaotic 
quantization. 




FIG. 3: Complex Polyakov line values from 4-dimensional 
quantum Monte Carlo simulation (red dots) and from (4+1)- 
dimensional classical Hamiltonian equation of motion (green 
dots) at aE™ ag /T = S 4 /h. 



